A BCS-BEC crossover in the extended Falicov-Kimball model: 
Variational cluster approach 
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We study the spontaneous symmetry breaking of the excitonic insulator state induced by the 
Coulomb interaction U in the two-dimensional extended Falicov-Kimball model. Using the varia- 
tional cluster approximation (VGA) and Hartree-Fock approximation (HFA), we evaluate the order 
parameter, single-particle excitation gap, momentum distribution ftmctions, coherence length of ex- 
citons, and single-particle and anomalous excitation spectra, as a function of U at zero temperature. 
We find that in the weak-to-intermediate coupling regime, the Fermi surface plays an essential role 
and calculated results can be understood in close correspondence with the BCS theory, whereas in 
the strong-coupling regime, the Fermi surface plays no role and results are consistent with the picture 
of BEC. Moreover, we find that HFA works well both in the weak- and strong-coupling regime, and 
that the difference between the results of VGA and HFA mostly appears in the intermediate-coupling 
regime. The reason for this is discussed from a viewpoint of the self-energy. We thereby clarify the 
excitonic insulator state that typifies either a BGS condensate of electron-hole pairs (weak-coupling 
regime) or a Bose-Einstein condensate of preformed excitons (strong-coupling regime). 

PACS numbers: 71.10.Fd, 71.27.+a, 71.35.-y 



I. INTRODUCTION 

The realization of excitonic insulators (EI) on the prox- 
imity of the semimetal-semiconductor transition was sug- 
gested about half a century ago.^""^ Because of the weak 
screening of the Coulomb attraction between the elec- 
trons and holes due to the small number of carriers, the 
electrons and holes may spontaneously form bound states 
(excitons), giving rise to the EI state. As a candidate 
for EI, quasi-one-dimensional lT-TiSe2 has been studied 
both theoretically and experimentally i^i^ TaNi2Se5 has 
also been studied by angle-resolved photoemission spec- 
troscopy (ARPES) measurements. It was reported that 
the valence-band top is extremely flat, and the material 
can be a new candidate for an EI of bound pairs between 
Ni 3c?-Se 4p holes and Ta 5d electrons. ^'^ 

From the theoretical point of view, the Falicov-Kimball 
model'^ extended by including a finite valence bandwidth, 
i.e., extended Falicov-Kimball model (EFKM), has been 
extensively studied in context of the EI or electric ferro- 
electricityi^iiS The EFKM contains the large bandwidth 
c-electrons (with a hopping integral tc and on-site en- 
ergy Cc), small bandwidth /-electrons (with a hopping 
integral tf and on-site energy e/), and a Coulomb inter- 
action (U) between c- and /-electrons. The ground state 
phase diagram of the EFKM in the weak-to-intermediate- 
coupling regime was obtained by the constrained path 
Monte Carlo (CPMC)<^ In the strong-coupling regime, 
the EFKM can be mapped onto the spin-1/2 Ising-like 
XXZ model with a uniform magnetic field. In that case, 
the spontaneous EI ordering corresponds to the spon- 
taneous magnetization in the XY-plane and its phase 
diagram was also determined.— The phase diagram of 
the EFKM is composed of three phases: the charge- 
density-wave (CDW) with staggered orbital order (SOO) 



phase, excitonic insulator (EI) phase, and band insula- 
tor (BI) phase. The CDW phase is characterized by 
the periodic modulation of the total density of c- and 
/-electrons, while the SOO phase is characterized by 
the periodic modulation in the difference between the 
c- and /-electron densities. The instability toward the 
CDW and SOO phases was studied in detail by Zenker 
et ali^ The EI phase is characterized by the sponta- 
neous c-f hybridization. The BI phase is characterized 
by the completely filled c- or /-band. Interestingly, the 
ground state phase diagram in the weak-to-intermediate- 
coupling regime obtained by a Hartree-Fock approxi- 
mation (HFA)^^ agrees quite well with that obtained 
by CPMC. On the other hand, excitation properties of 
EFKM are still of great interest. Finite-temperature 
phase diagram and electron-hole bound state formation 
in EFKM were studied by HFA.^^ Projector-based re- 
mormalization method (PRM) calculation on the one- 
dimensional EFKMii reported that incoherent parts of 
the single-particle excitation spectra, which are related to 
the dissociation of the excitons, appears especially in the 
BEC regime. Detailed studies on the dynamical excitonic 
susceptibility at finite-temperature calculated by use of 
the PRM^'^ and slave boson (SB) techniquc"'^^ confirmed 
that tightly bound excitons exist even above the criti- 
cal temperature for exciton condensation. The results 
strongly support the so-called excitonic halo suggested 
by Bronold and Fehske,^- where tightly bound excitons 
exist without condensation, and the scenario of the Bose- 
Einstein condensation of preformed excitons in the semi- 
conductor side. Thus the effects of electron correlations 
on the static and dynamic properties of this model are 
worth studying. 

In this paper, we study the EI state of the EFKM de- 
fined on the two-dimensional square lattice as a function 
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of the Coulomb interaction strength U . We employ the 
variational cluster approximation (VCA)il based on the 
self-energy functional theory (SFT)^^ at zero tempera- 
ture. The cluster perturbation theory (CPT)^^ is used 
to calculate the single-particle Green's functions. We 
also employ HFA to clarify the effects of electron cor- 
relation that can be taken into account in VGA. As far 
as we know, VGA has not been applied to the study 
of the EI state of EFKM. The advantage of VGA com- 
pared to HFA is that VGA can fully take into account 
static and dynamic effects of electron correlations within 
the range of a finite size cluster. So far, VGA and 
GPT applied to a variety of strongly correlated elec- 
tron systems, such as the half-filled Hubbard model with 
competing magnetic orders in two-dimensionaP*^ and 
three-dimensional systems^^ periodic Anderson model 
with the competition between magnetic ordering induced 
by the Ruderman-Kittel-Kasuya-Yoshida interaction and 
nonmagnetic Kondo screening^^^ multi-orbital system 
with spin-orbit coupling, XY-plane magnetic ordering;^ 
etc, and it turns out that the method is useful to discuss 
correlation effects on the symmetry breaking or single- 
particle excitation spectra, especially in the insulating 
state. Thus it is worth studying the EI state of the 
EFKM by applying VGA and GPT to investigate the ef- 
fects of short range correlations on the symmetry break- 
ing or single-particle excitations. 

We will first discuss the U dependence of the calcu- 
lated EI order parameter and single-particle gap. In the 
weak-coupling regime, as expected from the BGS theory, 
the single-particle gap is scaled well by the order param- 
eter, whereas in the strong-coupling regime, the order 
parameter rapidly decreases with increasing U . Then we 
will show the calculated momentum distribution func- 
tions as functions of U . In the weak-coupling regime, 
the momentum distribution functions behave like those 
in the BGS theory, whereas in the strong-coupling regime 
the momentum dependence of the momentum distribu- 
tion functions becomes weak, and the behavior is quite 
different from BGS theory. The coherence length of the 
exciton shows a shallow minimum at the crossover regime 
as a function of U . We further calculate the single- 
particle spectra, anomalous Green's functions, and den- 
sity of states in order to investigate the electron correla- 
tion effects on the quasi-particles. Thus our study will 
shed light on the BGS-BEG crossover^"' in the EI state. 

This paper is organized as follows. In Sec. [Ill we intro- 
duce our model and method of calculation. In Sec. IIIIl 
we present our results for the EI order parameter, single- 
particle gap, single-particle Green's function, anomalous 
Green's function, momentum distribution functions, and 
coherence length as functions of the Goulomb interaction 
strength U . Discussion on the efficiency of the HFA on 
this model and experimental implications are given in 
Sec. IIVI We summarize our work in Sec. |Vl 



II. MODEL AND METHOD 

A. Extended Falicov-Kimball model 

The Hamiltonian of the EFKM reads 
n = -ic ^(cfcj +H.C.) 4- (ec - 

-i/E(///^+H-c-) + (e/-M)E^»/ 

+U^nicnif (1) 

i 

where Ci (c|) denotes the annihilation (creation) operator 
of an electron on the c-orbital at site i and n^c = c|ci. 
tc is the hopping integral between neighboring sites of 
the two-dimensional square lattice and Cc is the on-site 
energy level of the c-orbitals. These are the same for the 
/-orbitals. U is the inter-orbital Goulomb repulsion be- 
tween electrons. The chemical potential /x is determined 
so as to maintain the average particle density n at half 
filling n — 1. Throughout the paper, we set h — ks = ^ 
and lattice constant a — I. We use = 1 as the unit 
of energy and we focus on the band parameter values 
Cc — 0,€f = —I, and tf ~ —0.3. 

The non-interacting tight-binding band structure and 
corresponding Fermi surface at half-filling are shown in 
Fig. [TJ At this parameter set, c- and /-band are over- 
lapping each other. The level difference e/ — Cc = — 1 
causes an imbalance between the c- and /-electrons den- 
sity, i.e., ric = 0.34 and nf — 0.66, where ric and nf 
are the average density of c- and /-electrons respectively, 
and we can see that Q (tt, tt) is not a nesting vector of 
the Fermi surface. Thus, in this paper, we do not con- 
sider any periodic modulations characterized by Q, such 
as the GDW phase, which is realized with small energy 
difference |ec — e/I and ric = rif = O.SJ^ 

B. Variational cluster approximation 

In order to analyze the EI state of EFKM, we apply 
the VGA.-'^^ Here we briefiy review the formulation of the 
VGA in order to make our paper self-contained. Follow- 
ing Potthofffi2. the grand potential functional is given as 

= J-p] -Trln(-G(7i + S) (2) 

where is the Legendre transform of the Luttinger- 

Ward functional ^[G]r^ Gq is the non- interacting 
Green's function, and we call S the trial self-energy. Tr 
represents the sum over fermionic Matsubara frequencies 
with temperature T and trace over the single-particle ba- 
sis. The explicit definition of Tr will be given later (in 
Eq. ([in])). The stationarity condition SHI'S] /Si: = 
gives the Dyson equation, and the functional gives the 
grand potential of the system at the stationary point 
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(0,0) 



FIG. 1: (Color online) (Top) The non-interacting tight- 
binding band structure of the c-orbital (dashed line) and /- 
orbital (dash-dotted line). Parameter values are tc = l,ec — 
0,6/ = —1, and tf — —0.3. The horizontal solid line rep- 
resents the chemical potential at half-filling. (Bottom) The 
non-interacting Fermi surface at half-filling. The momentum 
path (0, 0) ^ (tt, tt) — ^ (tt, 0) — s> (0, 0) is also shown (thick 
straight line). 



SFT— provides a way to compute by using the fact 
that the functional form of J-[S] depends only on the in- 
teraction terms of the Hamiltonian. Here we introduce a 
so-called reference system, which consists of disconnected 
finite-size clusters forming a super-lattice. Note that, be- 
cause the interaction term of the Hamiltonian Eq. ([T]) is 
local, the original and reference systems have the same 
interaction term. Therefore the functional form of 
is unchanged. The exact grand potential of the refer- 
ence system is given as fi' = T[E'] — Trln(— G'q ^ + S'), 
where S' is the exact self-energy of the reference system. 
Then, by restricting the trial S to S', we can omit the 
functional and obtain 



n[E'] = n' -Tiln{I -VG') 



(3) 



— Gg ^ represents 



where I is the unit matrix, V 
the difference of the one-body terms between the original 
and reference systems, and G' = (G'g ^ — S')~^ is the 
exact Green's function of the reference system. The size 
of these matrices are 2Lc x 2Lc, where Lc is the number 
of sites within a disconnected finite-size cluster. 



The trial self-energy for the variational method is gen- 
erated from the exact self-energy (or the exact Green's 
function) of the reference system. The Hamiltonian of 
the reference system is defined as 



— 7^ + Hpair + "Hlocalj 



H 



pair 



= A'^(cl/. + H 



^locai = jriic + Uif) , 



(4) 
(5) 

(6) 



where the Weiss field for the on-site electron-hole pair- 
ing A' and the orbital-independent potential e' are varia- 
tional parameters, which are optimized based on the vari- 
ational principle, i.e., {dn/dA',dn/de') ^ (0,0). Note 
that the solution with A' 7^ corresponds to the sponta- 
neous EI state, e' is introduced in order to calculate the 
average particle density n correctlyi^ Then we solve the 
ground-state eigenvalue problem 'H'lV'o) = £'o|V'o) of a 
finite-size cluster and calculate the trial single-particle 
Green's function by the Lanczos exact-diagonalization 
method. The Green's function matrix in Eq. ^ is de- 
fined as 



G'( 



G'^'iu) G'"{u) 



(7) 



where G"^'s are the Lc x Lc matrices. Each matrix ele- 
ment is defined as 



1 



LJ-W 
1 



n'-Eo 



(8) 



and these are calculated by the standard Lanczos tech- 
nique. The matrix V is given as 



V(K) 



T'=(K)-e'/ -A'/ 
-A'7 Tf{K)-e'I 



(9) 



where T"(K) is the inter-cluster hopping matrix for a- 
electrons. The matrix elements are given as T^CK) = 

taJ2yi.x^^^'^^i+xd^R+^,R'' where x denotes the neigh- 
boring sites of the z-th site and X denotes the neighboring 
clusters of the R-th cluster. 

Now all the physical quantities are diagonalized for 
Matsubara frequencies and super-lattice wave vectors, 
but not for orbitals and sites within a cluster. Thus Tr 
for a quantity A is written explicitly as 

TrA = T^e'"""+^ ^ J^A^{K,icv,,). (10) 

uin K a=c,f 1=1 

The K-summation is done in the reduced Brillouin zone 
of the superlattice. For numerical calculations of fl, the 
Matsubara-frequency sum is transformed to a contour in- 
tegral with complex Fermi function /(w) = l/(e"/"'" -I- 1) 
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by the theorem of residuum. Then the contour is de- 
formed to a path enclosing the real axis by use of the 
convergence factor. Finally we obtain an expression for 
the functional, 

n = n' - (( — Y Indct (7 - F(K)G'(cj)) . (11) 

'' K 

The single-particle Green's functions are calculated by 
CPT^^ with the optimized variational parameters. The 
CPT Green's function is defined as 



(7"^(k,..) = ^ 1^ Gl 



CPT,ij 



(k,cj)e" 



^ik(ri 



(12) 



where is the position of the z-th site within a discon- 
nected finite-size cluster and GcPT(k,t<j) = G'{uS)(l — 
V(k)G'(a;))~"'^. The wave-vector k can take arbitrary 
values in the Brillouin-zone. Here we define t/'^'^(k, w), 
Q^^ uj), and ^'^■''(k, iv) as the single particle c-electron, 
/-electron, and anomalous Green's functions, respec- 
tively. 

A cluster of the size — 8 (16-orbital) is used as a ref- 
erence system, thus the effects of statical and dynamical 
electron correlation within the cluster size are taken into 
account. Details of VGA can be found in Refs. [27ll28j . 



C. Hartree-Fock approximation 

It was reported that the ground-state phase diagram 
of the two-dimensional EFKM obtained by HFA quan- 
titatively agrees with that by GPMG^ in the weak-to- 
intermediate coupling regime."':^ To compare VGA results 
with HFA results, we briefly review the mean-field the- 
ory for EI state of this model.— Applying the HFA to the 
interaction term in the original Hamiltonian Eq.([T]), i.e., 
4cjlf. ^ {cU)fUMflh)cU-{f}c,)c\h-{c\h)flc., 
where (• • ■) denotes the ground-state expectation value, 
and switching to momentum space, we obtain the mean- 
field Hamiltonian 



-A 



Ck 

/k 



-A ek/ 

Eka = 2t ct (cos kx + cos ky) + ta " M + UUa , 
k 

u 



(13) 



where L is the number of lattice sites and Ck (/k) is the 
Fourier transform of Ci (fi). a = c, f represents the or- 
bital index and a denotes the other orbital of a, i.e., 
c = f and vice versa. The order parameter A describes 
the coherent exciton formation between c-electrons and 
/-holes. Here we assumed A is real without loss of gen- 
erality. Introducing the fermionic quasi-particles defined 



7k+ 
7k 



Uk \ I Ck 

Vk ~Uk / I /k 



(14) 



with + — 1 and diagonalizing the matrix in Eq. 
for each k, we obtain the mean-field Hamiltonian 

Hhfa = fe7k+^k+ + i?k 7k-Sk-) , (15) 



with the quasi-particle dispersion 
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^k =^(^kc + £k/)±Ve^ + A2, (16) 



■^k - 2(<^kc - fk/)- 



(17) 



Self-consistency equations for the particle density and the 
order parameter are 



k 

^E(-k/(i?k+) + -k 

k 

U V- 



(18) 
(19) 



A = -YukVk{l{E+)-f{E^)), (20) 

k 

respectively, where the quasi-particle density (7k ^7k) is 
replaced by the Fermi function f{E^) = l/(e-^k/^ + l). 
The coefficients are given as 



^ 1 



6 



v/eJTA^, 



Ck 



A 



^ -L A2' 



(21) 

(22) 
(23) 



The parameters ric, rif, and A are determined by solving 
the above equations self-consistently. 

Inverting the matrix {u — Hk), where Hk is the 2x2 
matrix in the mean-field Hamiltonian (|13p . we obtain the 
c-orbital, /-orbital, and anomalous Green's function 

e^rFA(k, uj) = , (24) 



UJ - ek, 



1 



ui-ekf 
a: 

^k/ LU-Hkc 



A2 ' 



(25) 



(26) 



{u! ~ ekc)(t^ - Ck/) - A2 ' 

respectively. From the imaginary part of each of the 
Green's functions, we obtain the c-orbital, /-orbital, and 
anomalous spectral function 

^W(k,^) = ul5{uj ~ E+) + vlS{u - E^ ), (27) 

4FA(k,^) = vlSiuj - E+) + ulSiiu - E^), (28) 

FHFA(k, uj) = UkVk {S{lu - E+) - S{lj - E^)) ,(29) 
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FIG. 2: (Color online) U dependence of the order parame- 
ter 2|A| and the single-particle gap Agap calculated by VGA 
(upper panel). HFA results for 2|A| and Agap are also shown 
(lower panel). The dash-dotted straight line indicates the 
single-particle gap in the normal (A = 0) state. 



respectively. 



We can see from the results that there are not only a 
lower bound of the Coulomb interaction strength Ud but 
also an upper bound Uc2 for the EI state. The obtained 
values are (C/d, Uc2) = (0.65, 6.6) for VGA and (C/ci, f/c2) 
— (0.66, 6.95) for HFA, respectively. The existence of 
the upper bound Uc2 seems to contradict to the case of 
the attractive Hubbard model, which has no Uc2'^^^ 
What happened aXU = Uc2 is that, the Hartree potential 
makes the /-band fully occupied and c-band empty, so 
that there is no Coulomb interaction between c-electrons 
and /-holes. Thus the system is simply a band insulator 
above Uc2- 

Note that, which band becomes empty or full at large 
U is determined by the particle density of each orbital 
at J7 = 0. In our case, ric < nf aX U=Q. Therefore the 
Hartree potential for c-band is larger than that for /- 
band. Thus, with increasing U , the c-band is pushed up 
rather than / band is, and finally c-band becomes empty. 

In the weak-to-intermediate coupling regime {U < 5), 
both the order parameter and single-particle gap increase 
with increasing U with the relation 2|A| ~ Agap (up- 
per panel). This result is consistent with the relation 
2|A| = Agap from HFA (lower panel). In the strong- 
coupling regime (U > 5), the order parameter rapidly 
decreases with increasing U but the single-particle gap 
remains open. If we can assume that the energy scale 
of the single-particle gap Agap and order parameter 2| A| 
may correspond to that of the characteristic temperature 
for the exciton formation (Tex) and critical temperature 
for the condensation of excitons (Tc), respectively, then 
the two temperatures should be comparable (Tex — Tc) 
in the weak-coupling (ECS) regime but may be well sep- 
arated (Tex S> Tc) in the strong-coupHng (BEG) regime. 
The BGS-BEG crossover may then be expected in this 
model although our calculations are done at zero tem- 
perature. 



III. RESULTS OF CALCULATION 

A. Order parameter and single-particle gap 

We first calculate the U dependence of the order pa- 
rameter for exciton condensation 

2A=^^(cl/. + H.c.) (30) 

i 

and the single-particle excitation gap defined as 

Agap = A*^ - (31) 

where /j.'^^^^ is the upper (lower) bound of the chemical 
potential. Calculated results by VGA (HFA) are shown 
in the upper (lower) panel of the Fig. [2l The factor 2 for 
the order parameter is introduced in order to compare 
with the single-particle gap, by analogy with the BGS 
mean-field theory. 



B. Momentum distribution function 

We then consider the c-electron, /-electron, and 
anomalous momentum distribution functions defined as 

n„(k) = / ^e""(k,z), (32) 

Fik) = / ^g^f{Kzl (33) 

respectively, where the contour integral path C< encloses 
the poles of the integrand on the real axis below the 
chemical potential. The results are shown in Fig. [31 Here 
we define, by analogy with the BGS mean-field theory, 
the "Fermi momentum" kp as 

nc(kF) = n/(kF) = 0.5. (34) 

Actually in HFA, this definition means that the single 
particle gap is identical to 2|A| and 2\F{'k)\ = 1 at kp 
(see Eq. ^ and Eq. ^ respectively). At [/ = 2, 
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F(k) shows the sharp peak at kp, indicating the exis- 
tence of weakly bound electron-hole pairs. At U — 6.5, 
k-dependence of the peak intensity of the anomalous 
Green's function is weak and F(\i) is spread out in mo- 
mentum space. Thus, in real space, small electron-hole 
pairs exist in the strong-coupling regime. 




-1.0 



f*™,/ y — — U=6.5 



V 



(0,0) 



(7I,7l) 



{n,0) 



(0,0) 



FIG. 3: (Color online) Calculated momentum distribution 
functions nc(k), n/(k), and F(k) a,t U = 2, U = 5, and 
U = 6.5. 



To see the U dependence of the momentum distribu- 
tion functions in more detail, we show the results along 
the (0, 0) (tt, tt) line in Fig. H) We can see that, in the 
weak-coupling regime, ric(k) drops sharply across kp and 
F(k) is peaked at kp. With increasing U, kp approaches 
(0, 0) because the Hartree potential for the c-electron re- 
duces the c-electron density, and -F(k) becomes broad in 
momentum space, indicating that the radius of electron- 
hole pairs becomes small in real space. When U reaches 
the crossover regime {U ^ 5), we have no kp and |i^(k)| 
decreases for all momenta with increasing U. This be- 
havior is consistent with the rapid decrease of |A| in the 
strong-coupling regime (see Fig. [2]). 



C. Coherence length 



We also evaluated the coherence length defined as 



^coh 



^ EklVkf(k)F 



(35) 



in order to see the spatial coherence of the excitons di- 
rectly. The k-summations were done with 100 x 100 k- 
points in the 1*** Brillouin zone. For VGA calculations, 
the derivative with respect to kx,y was evaluated by the 
4-point finite difference, while for HFA calculations, the 
analytical expression for Vki^(k) was used. Galculated 
results are shown in Fig. [5] In the weak-coupling regime, 
''coh is spread widely, about several lattice-constants and 
rapidly decreases with increasing U. Note that the calcu- 
lated result of Tcoh by VGA is considerably smaller than 




(n,n) 



U=l - 


- U^5 


U=2 - 


- U=6 


U=3 - 


- U=6.5 


(7=4 




(7l,7t) 



FIG. 4: (Color online) The momentum distribution functions 
nc(k) and 2_F(k) along (0, 0) — >■ (tt, tt) for various U calculated 
by VCA (upper panel) and HFA (lower panel). The /-electron 
momentum distribution functions n/(k) = 1 — nc(k) are not 
shown. 



that by HFA in the weak-coupling regime, especially for 
Tcoh >1. AtC/ = 2^3, Tcoh is already the size of 
the lattice-constant. A similar rapid decrease of the co- 
herence length of the Gooper pairs as a function of the 
Goulomb interaction strength was also reported in a de- 
tailed exact-diagonalization study on the attractive Hub- 
bard model. '^^ Furthermore, we find that rcoh has a shal- 
low minimum at C/ ~ 5 where the system is in crossover 
regime. This is because the denominator |_F(k)p in 
Eg. ([55)1 is largest in the crossover regime (see Fig. U). 
Then rcoh slightly increases again with increasing U. 



D. Single- particle spectra and density of state 

We also calculated the single-particle and anomalous 
excitation spectra 



A(k,w) = -^3 y g""(k,cj + i77) (36) 

(37) 



F(k,w) = 3e=-^(k,w + i77) 

TT 
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FIG. 5: (Color online) U dependence of the coherence length 
calculated by VGA (circles) and HFA (solid Hne). 



and the density of states (DOSs) 

P"H = ^E^"(k,a.) (38) 

k 

at [/ = 2 (BCS regime), U = 5 (Crossover regime), 
and U = 6.5 (BEC regime). The results are shown 
in Fig. m and Fig. [71 respectively. The HFA quasi- 
particle dispersion is also shown in the lower panel 
of Fig. [S] Note that, because of the artificial supercell 
structure introduced by the VGA, the spectra show ar- 
tificial gaps due to Brillouin zone folding. At U = 2, 
A(k, (jj) shows a small gap at the Fermi momentum kp 
defined in Eq. p4l) . F(k, w) shows sharp peak near kp 
and its intensity rapidly decreases as the momentum goes 
away from kp or the frequency goes away from the Fermi 
level fj,. At U = 5, the incoherent continua appear in the 
spectral function. We can see from the DOS and anoma- 
lous Green's functions that both the single particle gap 
and hybridization are large. At U = 6.5, A(k, ui) shows a 
semiconductor-like dispersion mainly due to the Hartree 
potential. The momentum dependence of the intensity 
of F(k, cj) is weaker than that a.t U — 2 and 5. Note 
that, although U is large, the incoherent part of both 
A(k, oj) and F(k, oj) is smaller than that at [/ = 5. The 
dispersion relation is well described by the HFA quasi- 
particle dispersion both in the weak- and strong- 
coupling regime. The reason will be discussed from the 
view point of the self-energy at Sec. |lVl 

IV. DISCUSSIONS 

A. Why HFA works well 

Here we consider, from the view point of the self- 
energy, why the HFA is successful not only in the weak- 



coupling regime but also in the strong coupling regime for 
EFKM. For simplicity, we neglect the order parameter A 
and Weiss field for electron- hole pairs A'. 

Using the spectral representation^^ the self-energy can 
be written as 

S(k,..) = 5k + E^T^' (39) 

where is the Hartree potential^ Ck,iy is the v-th pole 
of the self-energy, and (Jk.i/ is the corresponding spectral 
weight. In HFA, the Hartree potential is taken into ac- 
count but the second term (frequency dependence of the 
self-energy) is neglected. Note that EFKM defined in 
Eq. ([T]) is nothing but the asymmetric Hubbard model. 
Therefore we can apply the sum-rule for the self-energy of 
the Hubbard model (see Ref. fsH and Appendix). Then 
total weight of the self-energy 'neglected' in the HFA is 

1 f°° 

lim / da;3I](k, a; + ir?) = > a-k v = U^ricnt . 

7r')->o+y_^ > ^ ■ 

(40) 

Here the half- filling condition Uc + Uf — 1 is used. The 
U dependence of ric, nj, and U^ncUf calculated by the 
HFA is shown in Fig. El We can see from the result that, 
with increasing U, the Hartree potential causes particle 
number imbalance of the c- and /-orbitals, and makes 
ncUf smaller. Thus U'^ricnf has a maximum and decrease 
again with increasing U. 

Moreover, we can see from the HFA Green's functions 
Eq. ([M]) that, if the order parameter A is finite, the 
single-particle gap can be generated from the hybridiza- 
tion gap in the weak-coupling regime. Thus HFA works 
well on the EFKM both in the weak- and strong-coupling 
regime. 

The crucial differences between VGA and HFA will ap- 
pear in one- and two-dimensional system at finite tem- 
perature, where the spontaneous symmetry breaking is 
absentj^^i^ Actually, it was reported that the critical 
temperature for exciton condensation of the EFKM eval- 
uated by the SB technique is lower than that by the 
mean-field due to the effects of electron correlations)^ 
Moreover, the effects of the spatial fluctuations, which 
are also completely neglected in the mean-field theory, 
will also tend to destroy the ordering. 

B. Experimental implications 

Since the EI order parameter is not necessarily identi- 
cal with the single-particle gap, especially in BEG regime, 
an experimental evidence for the realization of EI should 
be signaled as the spontaneous hybridization between the 
valence and conduction bands. Generally, ARPES exper- 
iments observe the imaginary part of the single-particle 
Green's function filtered by the dipole matrix element 
and the Fermi functioni^ The matrix element effects are 
determined by the selection rule from symmetries, pho- 
ton energy dependence of the cross section, etc. Gon- 
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FIG. 6: (Color online) Calculated anomalous spectra -F(k, uj) (upper panel) and single-particle excitation spectra ^(k, oj) (lower 
panel) aX U = 2 (left), U = 5 (center), and U = 6.5 (right). In the lower panel, the HFA quasi-particle dispersion is also 
shown (solid line). The artificial Lorenzian broadening r; = O.I is used. 



sideration about the matrix element effects becomes im- 
portant for observation of Zhang-Rice singlet states in 
cuprate by ARPES.''^ If one can resolve dominant or- 
bitals for the valence and conduction band by use of ma- 
trix element effects, the spontaneous hybridization be- 
tween the valence and conduction bands can be observed 
by ARPES as the difference in the spectral intensity be- 
tween above and below the critical temperature. That 
is, the spectral intensity from the dominant conduction- 
band orbitals will be transferred to the valence-band top 
below the critical temperature. Therefore, experimental 
analyses of temperature dependence of the hybridization 
between the valence and conduction bands are desired. 

The recent ARPES measurements on quasi-one- 
dimensional Ta2NiSe5 were done with the photon ener- 
gies hf = 10 eV and hiy = 23 eVi^ The experimental 
results of the energy distribution curve (EDC) at T = 40 
K showed that the EDC intensity is large near F with 
hiy = 23 eV, but small with hi/ — 10 eV. From the cross 
section table of atoms, "^^ Ta 5d weight should be large 
for hiy — 23 eV, while Se 4p weight should be large for 
hiy = 10 eV. The result implies that the spectral weight 
of the conduction-band Ta 5d orbitals is transferred to 
the valence top due to the hybridization. The similar 
spectral-weight transfer can be seen in our calculated re- 
sult for the orbital-resolved excitation spectra shown in 



Fig. [HI We can see from the results that, near k = (0, 0), 
the spectral weight of the dominant conduction-band or- 
bital (c-orbital) is transfered below the Fermi energy. 
This spectral-weight transfer also can be seen in the cal- 
culated DOS at J7 = 5 (middle panel of Fig. [71). Thus, 
the temperature-dependent photoemission spectroscopy 
measurements with various photon energies are desired 
to identify the orbital character of the band structure 
and estimate the hybridization between the valence and 
conduction bands in Ta2NiSe5. 



V. SUMMARY 

In this paper we have analyzed the excitonic insula- 
tor (EI) state of the extended Falicov-Kimball model 
(EFKM) by using the variational cluster approximation 
(VGA) and Hartree-Fock approximation (HFA) at zero 
temperature. We have calculated the EI order param- 
eter, single-particle gap, momentum distribution func- 
tions, coherence length, and single-particle Green's func- 
tions, as a function of the Goulomb interaction strength 
U. 

In the weak-coupling regime, we found that the mag- 
nitude of the single-particle gap Agap is almost compa- 
rable to that of the order parameter 2|A|. This indicates 



co-iu 



U 




FIG. 7: (Color online) Calculated density of states a,t U — 2 
(top), U — 5 (middle), and U — 6.5 (bottom). The verti- 
cal line represents the Fermi energy. The artificial Lorenzian 
broadening rj = 0.1 is used. 



that the electron-hole pair formation and its condensa- 
tion may occur simultaneously, like Cooper pair forma- 
tion and its condensation in BCS theory. The quasi- 
particle dispersion obtained by the single-particle exci- 
tation spectra ^(k, w) is well described by mean-field 
theory. The Fermi momentum kp is defined from the 
momentum distribution function. The anomalous excita- 
tion spectra i^(k, cj) showed that its spectral weight dis- 



FIG. 8: (Color online) U dependence of ric (dashed line), 
Uf (dash-dotted line), and U'^ricnj (solid line) calculated by 
HFA. 
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FIG. 9: (Color online) Calculated single-particle excitation 
spectra A(k, a;), yl'^(k, w), and (from left to right) at 

(7 = 5. The horizontal line represents the Fermi energy. The 
artificial Lorenzian broadening 77 = 0.1 is used. 



tributed mainly near the Fermi energy /i. The anomalous 
momentum distribution function F{\<l) is peaked at kp. 
Reflecting this, the coherence length of the exciton spread 
wide for several hundred lattice spacing. This indicates 
that the system is in the BCS-like weakly-bound exci- 
ton condensation state. With increasing [/, the momen- 
tum dependence of the condensation amplitude F(k) be- 
comes weak, and the coherence length decreases rapidly. 
A(k, uj) show incoherent continua in their high-frequency 
part in the intermediate-coupling regime. In the strong- 
coupling regime, the energy scale of the order parameter 
and single-particle gap became separated. This result in- 
dicates that the binding energy of the electron- hole pairs 
(excitons) is larger than the energy scale of the critical 
temperature where exitons may obtain coherence. The 
Fermi momentum kp became ill-defined. Accordingly, 
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the Fermi surface plays no roles and F(k) is widely spread 
in momentum space, and the coherence length is smaller 
than the lattice constant, indicating that the system is 
in the BEC-like condensation state of strongly bound 
electron-hole pairs. Moreover, we found that HFA works 
well not only in the weak-coupling (small U) regime, but 
also in the strong-coupling (large U) regime. The reason 
was clarified from the view point of the self-energy. Fi- 
nally, we discussed the spectral feature of the EI state of 
the EFKM and gave experimental implications for pho- 
toemission spectroscopy measurements. 
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gives the Green's function with the form 

G'cr(k,a;) = (w - e^a - T,c,(k,uj)y 



(43) 



Then we substitute the spectral representation of the self- 



energy^ 



into (|43p and take the high-frequency expansion. 



(44) 



^ /, s _ 1 Cko- + 5ko- , <^ka,iy + (fkcr + ffkcr) 

(j"(T(,lt, — \ K 1 ^ 

+0 . (45) 

Substituting this expression into Eq. (|42|) and using the 
theorem of residuum, we obtain 

^^L = Xl^k<T,^. + (eka+5ko•)^ (46) 
The other expression for the second moment is given 



Appendix 

Here we derive the sum-rule for the self-energy used in 
Eq. (gOI). The EFKM or the asymmetric Hubbard model 
in momentum space is give as 



■H = ^ eko-c^Ck^ ^\Y^ 4+qtCkt4'-qiCk';. (41) 

kk'q 



kCT 



Here we consider the second-moment of the Green's 
function of the electron with spin u and momentum k, 
which is defined as^ 



ML = /^-^G.(k,..), 



(42) 



where Go- (k,a;) is the single-particle Green's function of 
the electron with spin cr, and the integral path encloses 
all singularities of the integrand. The Dyson equation 



AC = ({[ck.,H],[H,4j} 



(47) 



where {•■•}+ denotes the anticommutator. Calculating 
the (anti) commutators on the right-hand-side, we obtain 



(48) 



where a denotes the opposite spin direction of a. Now 
we use the fact that g\^„ is the Hartree potential)^ i.e., 
ffkcr = Una. Then comparing Eq. with Eq. (j^S]), we 
obtain 



^ CTkrr,!/ = U'^ng{l - rig) 



(49) 



The right hand side does not depend on the momentum k 
or dispersion eko-. By replacing a =t, 4- to c, / and using 
the half- filling condition ric + Uf = 1, we obtain Eq. (|40p . 
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